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1. Introduction 

The free parameters of QCD can be fixed by matching a corresponding number of hadronic 
observables at low energies. Then, predictions can in principle be made for any other observable 
both at low and high energies. A quantity of particular interest is the strong coupling constant a s . 
A commonly used reference value is the MS coupling of dimensional regularisation at the scale of 
the Z-boson mass. This poses a particular challenge for the lattice formulation, due to the large 
scale differences involved, which cannot be accommodated on presently affordable lattices. This 
difficulty can be overcome by using finite size techniques [JI| ||]. The essential ingredient is the 
non-perturbative definition of a running coupling, g (L), which runs with the linear extent L of the 
space-time volume. The Schrodinger functional provides a framework for such a definition, where 
g(L) is defined through the response of the system to a constant colour electric background field. 
See [^] for a detailed explanation and [f|] for the definition of the SF for QCD both in the continuum 
and on the lattice with Wilson quarks. 

Within this framework, the step scaling function, (SSF), plays a fundamental role. It can be 
regarded as a discretised version of the Callan and Symanzik j8 function, and can be used to study 
the evolution of g (L) with energy. Such studies have been carried out for different numbers of 
flavours [§, g— §|. We present here the step scaling function and the running coupling with four 
flavours of staggered quarks. In ^ an analogous work using 0(a) improved Wilson quarks was 
presented. The agreement of the results in the continuum limit can thus be regarded as a test of 
universality. 

This write-up is organised as follows. Section 2 reviews the basics of the SF on the lattice with 
staggered fermions, gives a definition for the renormalised running coupling, and revisits the basics 
on finite size techniques. In section 3, the 0(a) improvement for our setup is revisited. Section 4 
presents the details of the simulations and the data analysis. In section 5, we present our results 
and we finish with an outlook to future work. 

2. Schrodinger functional, coupling constant and finite size techniques 

The SF is a useful tool to study the scaling properties of QCD. Here we will introduce it 
briefly. See ^ |lO|-|T2|] for more details and [13-15] for the set up of the SF with staggered 



fermions. It can be regarded as the Euclidean time evolution kernel for going from a state at time 
xq = to another state at time xq = T . Using the transfer matrix formalism, it can be expressed as 
a path integral with fields satisfying periodic boundary conditions in space and Dirichlet boundary 
conditions in time. Homogeneous boundary conditions are imposed on the fermionic fields, 

P + Y\ Xo =o = = <^+L=r> P M X0 =T = = V^-L=o> (2-1) 
with P± = i(l ±)f)) and the spatial components of the gluon fields satisfy, 

M*)\ Xo=0 = Ck, A *WU=r = c *- ( 2 -2) 
The SF is then a functional of these boundary fields, 

3f[C',C\ = J ^{A,\i/,W SlA ' V,f] - (2-3) 

The choice of the boundary fields is largely arbitrary. Following we choose Abelian and spa- 
tially constant fields such that the absolute minimum of the action is unique up to gauge trans- 
formations. A judicious choice of the boundary fields ensures that the absolute minimum of the 
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action is unique up to gauge transformations and yields small lattice artefacts in the renormalised 
coupling. 

2.1 Definition of the coupling constant 

The SF allows us to define a renormalised coupling constant g 2 non-perturbatively, easily 
computable on the lattice and in perturbation theory and with reasonably small lattice artifacts. 
Since the induced background field is unique, it is possible to unambiguosly define the effective 
action of the SF, i.e., r[g] = _^ {c c% (2 4) 

The boundary fields are parametrised by the scale L and T7, v, two dimensionless real parameters 
A renormalised coupling can be defined through, 



r' 



dr 



T)=V=0 



8 2 (L) 



(2.5) 



g 2 (L) is going to run with the spatial box size L which plays the role of the scale. The normalisation 
constant is chosen so that the renormalised coupling coincides with the bare coupling at tree level 
in the perturbative expansion. 

2.2 SF with staggered fermions 

The SF with staggered fermions requires lattices where the time extent T/a is odd and the 
spatial extent L/a is even. [13, 14]. This is illustrated in figure [l] in a two dimensional sketch, 



as well as one possible reconstruction of the staggered fermions. See [16] for a more detailed 
description on the fermionic reconstruction. The continuum limit for g 2 is usually taken setting 
T = L already for finite values of a. In order to define the continuum limit in our case, we are 
obliged to adopt modified conventions, [O]. Lattices with T = L±a are interpreted as having 
physical time extent T' = T + sa with s = ±1, so that the condition T' = L provides us with 2 
regularisations to our problem. This modification forces us to have the 0(a) improvement revisited. 
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Figure 1: SF with staggered fermions, 2-dimensional sketch. Left, T = L + a right, T = L — a. Thin (green) lines 
represent the lattice where the one-component staggered fermions live, thick (orange) lines represent the effective lattice 
for the reconstructed fermions and the brown dots stand for the sites where the reconstructed fermions live. 

2.3 Finite size techniques 

We want to compute the scale evolution of the renormalised coupling. We introduce then the 
step scaling function in the continuum theory, 

a(u)=g 2 (2L)\ u=f(L) . (2.6) 

On the lattice, it can be obtained as the continuum extrapolation from a sequence of pairs of lattices 
with sizes L/a and 2L/a, 
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a(u) = lim£(w, a/L). (2.7) 

The procedure is repeated for a range of u values in [g 2 (L m i n ),g 2 (L mia )], 

w =g 2 (L min ), u k = a(uk-\) = g 2 (2*L min ), k = \,2... (2.8) 

After 7, 8 steps, energy differences of 0(100) are bridged. At sufficiently large energies, perturba- 
tion theory can be applied to relate the SF coupling g 2 (L) with a perturbatively defined coupling, 
e.g. g 2 At low energies, the connection with physical units can be established through the com- 
putation of a hadronic quantity, e.g., F n L msx . 

3. 0(a) improvement 

Due to computer power limitations it is advisable to construct an action that is 0(a) improved 
so that the dominating lattice artifacts are cancelled. Following Symanzik's improvement pro- 
gramme [ |T7| ] this can be done by adding irrelevant local counterterms to the action, monitored by 
adjustable coefficients. These coefficients admit an expansion in perturbation theory. In the SF 
framework additionally to the volume counterterms there exist boundary counterterms. 

We want our observable, F' to be 0(a) improved up to one loop in perturbation theory. The 
action is taken to be 

S[U,x,X] =s g [u]+s f [u,x,x], 

W = ^£w(p)tr{l-[/(p)}, (3.1) 

Sop 

A T ~" 1 

Sf[U,x,x] = a * E L [^Uti{x)x{x+a(i)-X^U^x-a{i)x(x-afi)\ , 

xo=ax,fl ^ a 

where the sum over the gauge fields runs over all the oriented plaquettes p and w(p) are weight fac- 
tors that take the value 1 except for the boundary plaquettes, where they take the values c t (go) for 
the time-like plaquettes attached to the boundaries and lc s (go) for the spatial plaquettes. Concern- 
ing the fermionic part, x,X represent the one-component fermionic fields, homogeneous boundary 
conditions have been assumed for the fermions, a constant phase factor A^ 1 has been included, and 
Tfy = (— i)£v</i-W fl ^ the usual phase factors of the staggered fermions. 

There are no 0(a) effects arising from the bulk, so that the staggered SF will be 0(a) improved 
by including a couple of boundary counterterms. The only pure gauge counterterm relevant in our 
context takes the form, 

c?(gg)tr{F«F«}, c?(g§) = c f (0) < +g 2 (c^ +N f c^) + O(g 2 ), (3.2) 

where the superscript s = ±1 stands for the two possible regularisations. The coefficients above 
have been calculated for our set up obtaining, 

(0), = 2 4°' 1)[ =0.0274(2), c, (M)l =0.0077856(4), 

2 + s' d°' 1U = -0.4636(6), 4 M) -' =-0.0266(8). 

Concerning the fermionic part, there is also a boundary counterterm which, in terms of the single- 
component fields, takes the form [Jl6[|, 



1 We have chosen Ao = 1 , = e' e * = e l 5 . It leads to a smaller condition number of the free fermon matrix ni]| . 
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[d s (g 2 ) - 1] [A fc C4(jc)^(x+a^)- A M Uj(x-ol)^(jc-ol)] } VQ=fl 

x,& 

Kfeo) - 1] LI^MZW [A^(x)*(x + ^) - A^(x-a£)x(x-a£)] }^ =r _ fl (3.4) 

x,& 



The tree level value of di°^ was found to be di°^ = 1 + 1 Ql8|]. 



4. Simulations and data analysis 

In order to carry out our computations, we have made use of a customised version of the code 



offered by the MILC collaboration [19p, where the 0(a) improvement presented in the previous 
section was implemented. The simulations have been run for lattice sizes L/a = 4,6,8, 12, 16 and 
s = ±1. The statistics range from 60,000 measurements to 160,000 (the quantity v, [Q] has also 
been measured). The data analysis has been performed by using Uwerr . m p0(|. 

Instead of tuning the values of j8 = 2N c /g^ for different L/a so that they correspond to the 



same value of the renormalised coupling, we have followed the procedure proposed in [21]. We 
measure g 2 {L) for a set of values j6,L/a and generate an interpolating function. This function is 
then used to tune j8. The interpolation function takes the form, 

1 /3 ' ( 2N\ il 



g 2 (P,L/a) 2N fi ' V JS 
The data have been fitted by making use of the least squares method. The interpolated data 
inherit two sources of errors, statistical and systematic. 

5. Results 

Since we have two regularisations at our disposal, it is possible to perform an analysis of the 
lattice artifacts of our data. A line of constant physics is defined by the coupling in one regulari- 
sation and we evaluate the coupling for the other regularisation at the same values of j8 . One- loop 
perturbative cutoff effects can be subtracted, by defining the quantity, 

m± = u± x [l +u T (mf — nif)] , (5.1) 

where mf is the coefficient in the perturbative expansion, g± = g$ + mfg^ and ± stands for s = ± 1 
In Figure || we show the lattice artifacts for u and u^ l \ We have also measured the continuum 
extrapolation of the step scaling function. If we choose one regularisation to fix the physics, the 
lattice step scaling function can be computed in the two available regularisations. The continuum 
limit has to be shared. Perturbative effects can be subtracted, by defining E, , 

EW( M ,a/L) = .^"'f^ , 8 = W}_£W = 8 lia / L ) u + S2 (a/L)u 2 + Otf). (5.2) 
1 + o\{ajL)u y > 

In figure ^ we can see the continuum extrapolation for some of our data. The fit has been done in 
the following way, 

Ei {u s ,a/L) = o m i xed (u s )+Ai(a/L) 2 , Z-\(u s ,a/L) = a mix ed(«s) +A 2 (a/L) 2 , (5.3) 

where s stands for the regularisation chosen to fix the physics. The fit was performed excluding the 
data for L/a = 4. An analogous fit was done for L±\(u s ,a/L). 

We interpolate the data of the continuum extrapolation of a(u). In Figure || we present such 
interpolation, plotting o{u)/u vs u, together with the perturbative approximations to 1,2 and 3 loop 
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Cont.ext. of g 2 (a/L* =-1) at fixed u = g 2 (a/L,l) 



Cont.ext. of g 2 {a/L,s =1) at fixed u = g 2 (a/L,-l) 
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Figure 2: Analysis of the lattice artifacts of our system. Diamond (blue) points are the values of g 2 (L/a). Dashed 
lines correspond to the fit to these data and asterisks (green), the continuum extrapolation. Circles, (magenta) are the 
values for the same data after performing the perturbative subtraction, and the dotted lines their fits. Squares (orange) 
represent their continuum limit (displaced from the origin). The solid horizontal (red) lines are the lines of constant 
physics, given also by a (red) triangle slightly displaced from the origin. L/a = 4 are not included. 



Continuum extrapolation of the ssf, s = 1 , 2p fitted 



Continuum extrapolation of the ssf, s = -f , 2p fitted 




Crossed, square 111 




G 0.01 0.02 0.03 0.04 0.05 0.0S 0.07 

(a/L) 2 



0.03 0.04 

(a/L) 2 



u=0 7118 
0.05 0.06 0.07 



Figure 3: Continuum limit extrapolation of the step scaling function. Diamond (blue) points are the values of the step 
scaling function using data from the same regularisation and the dashed lines represent the fits. Circles (red) represent 
the lattice step scaling function taking data from two regularisations and the solid lines their fits. Asterisks (magenta) 
are the continuum extrapolations <j(u) of the lattice step scaling function. The renormalisation prescriptions, (values of 
u) are explicitly given in the plots. The graph on the right correspond to the regularisation s = I and the one on the left 
to s = -1. 

in PT. Our fit function was a polynomial of degree 6 where the first coefficients where set to the 



perturbative coefficients up to 2 loop [ |12| | in PT. The values corresponding to the largest couplings 
were not included in the fits. In the plot, we show the results obtained from the two different 
regularisations. They are correlated and therefore they can not be simultaneously used in a fit. The 
differences between the two can be regarded as systematic errors. 

6. Conclusions and outlook 

We have computed the SSF of the QCD coupling in the SF scheme with 4 flavours of massless 
staggered quarks. Unfortunately, the discretisation errors are fairly large. Some more effort will be 
put into the data analysis, and the A parameter still needs to be computed. The results are in rough 
agreement with data obtained with Wilson quarks [Sh, but a detailed comparison is still needed. 
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Figure 4: Step scaling function a(u). The dotted-dashed (red), dashed (blue) and solid (magenta) lines represent 
the perturbative 1-loop, 2-loop and 3-loop a(u). Diamonds (red) represent the extrapolated a(u) from the s = 1 reg- 
ularisation and circles (green) from the s = —1. Their fits (largest value excluded) are given by the solid thick (black) 
lines. 
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